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In a class of frustrated magnets known as spin ice, magnetic monopoles emerge as classical 
defects and interact via the magnetic Coulomb law. With quantum-mechanical interactions, 
these magnetic charges are carried by fractionalised bosonic quasi-particles, spinons, which 
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can undergo Bose-Einstein condensation through a first-order transition via the Higgs mech- 
anism. Here, we report evidence of a Higgs transition from a magnetic Coulomb liquid to 
a ferromagnet in single-crystal Yb2Ti207. Polarised neutron-scattering experiments show 
that the diffuse [111] -rod scattering and pinch- point features which develop on cooling are 
suddenly suppressed below Tc ~ 0.21 K, where magnetic Bragg peaks and a full depolari- 
sation of the neutron spins are observed with thermal hysteresis, indicating a first-order fer- 
romagnetic transition. Our results are explained on the basis of a quantum spin-ice model, 
whose high-temperature phase is effectively described as a magnetic Coulomb liquid, while 
the ground state shows a nearly coUinear ferromagnetism with gapped spin excitations. 

Introduction 

Emergent magnetic monopole^have recently been identified in many condensed-matter systemd^. 
In particular, recent experiments^ have shown that the low-temperature state of dipolar spin- 
j^gUOinJ materials, Ho2Ti207 and Dy2Ti207, can be effectively described as a classical magnetic 
Coulomb phase, which hosts these emergent magnetic monopoles obeying the magnetic Coulomb 
lav^S Theoretically, this has been understood from a classical model where (111) Ising mo- 
ments, i.e., pseudospin-1/2, interact mainly through a magnetic dipolar interactionPI. This classi- 
cal Coulomb phase is characterised by dipolar spin correlations with a power-law deca)^^!!^. 

Quantum variants of spin iceP^S] j^^yg ^i^^ hccn found in the magnetic pyrochlore oxides 
RE2TM2O7 (ref. IT) with the rare-earth (RE) ions Yb (refs [I8H25]) and Pr (refs [I5l[l6l|26l|27]) 
and transition-metal (TM) ions such as Ti and Ir, or Sn. In this quantum version of spin ice, the 
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magnetic monopole charges are expected to be carried by fractionalised pseudospin-1/2 quantum- 
mechanical quasi-particles, spinons, that acquire kinetic energy by hopping via a pseudospin-flip 
exchange interaction, if they remain deconfined with an energy gapPS. Thus, these systems provide 
a unique playground to study the intriguing quantum dynamics of monopolar spinons. Quantum 
spin ice actually harbours the analogous "quantum electrodynamics", where the "electric field" 
and "electric charge" describe the Ising magnetic moment and the associated monopole charge, 
respectively'^. In particular, Yb2Ti207 provides a prototypical system described by an effective 
pseudospin-1/2 quantum spin-ice model^l^ showing strong transverse quantum fluctuations of 
magnetic dipoles as a result of the following three characteristics: First, the crystal field structure 
of the Yb^+ ion is characterised by a ground- state Kramers doublet well- separated from the first 
excited doubleP^. Second, the Curie-Weiss temperature is weakly ferromagnetic with 6cw ~ 
0.53 K (ref. |24)) as in dipolar spin ice'". Third, the (7-tensor has planar components, = 4.18, 
that are larger than the (111) Ising components = 1.77 (ref. 122)) . enhancing the contributions to 
the magnetic neutron-scattering cross-sections from planar fluctuations. 

An early specific-heat measurement on a powder sample of Yb2Ti207 revealed a sharp sin- 
gularity at 0.21 K, signalling a first-order phase transition^'. A previous neutron- scattering exper- 
iment, on a single crystal in which we have subsequently shown that the specific-heat anomaly is 
reproduced (Supplementary Figure SI), found an almost coUinear ferromagnetic structure with an 
ordered Yb^+ magnetic moment of l.l/is, a uniform magnetisation along [100] and a slow relax- 
ation of the spins over times of the order of two hours well below Tq ~ 0.21 K (ref. |241) . These 
results may indicate the formation of an energy gap in the spin excitations below a first-order phase 
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transition to the ferromagnetic state. We have shown from specific-heat measurements (Supple- 
mentary Figure SI) and extended X-ray absorption fine structure (EXAFS) experiments (Supple- 
mentary Figure S2) on three different single crystals (see Methods), that the long-range ferromag- 
netic order in Yb2Ti207 is quite sensitive to Yb deficiency. In particular, our analysis indicates that 
as the Yb deficiency increases, the specific-heat anomaly becomes significantly broader and the 
sharp anomaly is lost (Supplementary Figure SI). This provides an understanding of the experi- 
mental results reported for different single-crystal and poly crystalline samples of Yb2Ti207, which 
suggest magnetic correlations remain short-range and dynamic even below 0.21 K (refs [T8ll23ll25l) . 
The variation in the experimental findings on different samples suggests the proximity of this mate- 
rial to a quantum phase transition between a quantum magnetic Coulomb phase (U(l) spin liquid 
and a magnetically ordered phase (Fig. 1 a). The magnetic ordering in this case can be under- 
stood as a Bose-Einstein condensation of monopolar spinons within the Higgs mechanisnP^*^, 
which gaps out all the soft spin excitations and takes a first-order charactei'^^^, as in the case of 
normal-to-superconductor transition^. This first-order quantum phase transition should extend 
to finite temperatures^^. Suppose that the transition temperature is much lower than an energy 2 J 
required for thermally and dynamically exciting monopolar spinons and that the high-temperature 
phase does not suffer from a static and thus non-uniform proliferation of monopole charges that 
causes a spontaneous symmetry breaking. Then, the high-temperature phase may be described as 
a magnetic Coulomb liquid at temperatures lower than 2 J, as in dipolar spin ice. Furthermore, 
the low-temperature phase transition can be effectively described as a Higgs transition of magnetic 
monopoles, like that of electric monopoles in the classical caseP'. The minimal requirement for 
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the identification of this Higgs transition of magnetic monopoles is the observation of a first-order 
phase transition separating a high-temperature Coulomb phase that may be characterised by pinch- 
point featureP'^'^'^ and a low-temperature magnetically ordered phase that exhibits finite planar 
components of pseudospins and gapped spin excitations. 

In this letter we report evidence of a Higgs transition from a magnetic Coulomb liquid to a 
ferromagnet, from data which were collected using low-temperature polarised neutron-scattering 
experiments on a single crystal of Yb2Ti207. We observe that well below Tq ^ 0.21 K the incident 
polarised neutron beam is completely depolarised. This is ascribed to the formation of macroscopic 
ferromagnetic domains. The first-order nature of the phase transition is revealed via the observa- 
tion of a thermal hysteresis around Tq. A nearly coUinear ferromagnetic structure consistent with 
previous neutron-scattering experiments on the same single crystai^is identified theoretically us- 
ing an effective pseudospin-1/2 quantum spin-ice model on a microscopic basi^^^"^^"^. The model 
also reproduces our polarisation-dependent magnetic neutron-scattering data above Tq, showing 
dipolar spin correlations with a remnant of a pinch-point singularity as a signature of an unusual 
finite-temperature paramagnetic phase which in practice can be described as a classical magnetic 
Coulomb liquicPH™. 

Results 

Diffuse scattering and pinch-point features above Tc Figure 2 a shows the experimental mag- 
netic diffuse neutron-scattering cross-sections collected at 0.3 K in the energy-integrated diffrac- 
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tion mode with an integration over the energy transfer up to 3.64 meV. The background has been 
subtracted (Supplementary Methods). Here, the incident neutron spins were polarised along the 
Z direction, perpendicular to the [hhV) plane, allowing us to probe the scattering intensities in the 
spin-flip (SF) and the non-spin-flip (NSF) channels as well as for the sum of these two channels. A 
broad intensity from (000) through (111) to (222) is identified in the sum (Fig. 2 a, right), as pre- 
viously reported in unpolarised neutron- scattering experiments'*^. This feature is dominant in the 
NSF channel (Fig. 2 a, middle). The [lll]-rod scattering is accompanied by a branch from (111) 
to (220) in the SF channel (Fig. 2 a, left), and by another from (222) towards (004) in the NSF 
channel (Fig. 2 a, middle). The anisotropic nature of the spin correlations is clearly seen around 
(111) in the SF channel (Fig. 2 a, left), suggesting a remnant of the pinch-point singularity^ that 
has been observed as a sharper singularity in dipolar spin icel^^. In order to examine this structure 
in more detail, we have performed similar measurements with a higher resolution at 0.3 and 1 .5 K. 
Figures 2 d and e show the neutron-scattering cross-sections for the SF channel along [11 — 1] and 
[111] cuts, respectively, through the (111) point. Whereas the intensity along [111] is featureless 
and simply increases in magnitude on cooling, the intensity along [11 — 1] clearly shows a growth of 
the correlation length. This provides direct evidence of the remnant of the pinch-point singularity. 

These polarisation-dependent anisotropic neutron- scattering results are reproduced theoret- 
ically as follows. The Yb'^+ 4/ magnetic moment at a site r is described with the ionic Kramers 
doubletPand thus the pseudospin-1/2 operator 5,. = (S'^, S*^, S^) as rh^ = {9i.{S%e% + S^e^) + g^S^e^^j 
with the locally defined coordinate frame (e^, e^, e^) (Supplementary Methods), where the Z axis 
points along the (111) direction. The effective Hamiltonian is comprised of the magnetic dipolar 
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interaction, 



{r,r') 



rfir ■ rhr' (rhr ■ Ar)(Ar ■ rh^/ 



(Z\r) (Z\r)5 
with Ar = \Ar\, and a nearest-neighbour anisotropic superexchange interactiorPEll^ 



(1) 



T n.n. 

Hse = -T. [SrSr' + '^'5+5,^ + gc^^^--' 5+5+ + Ke'^^^ {S',S+ + S^S',,)] + h.c, (2) 

with the bond-dependent phase factor (f)ry (Supplementary Methods). Here, J is the coupling con- 
stant for the Ising interaction found in spin-ice system^^SlllIl^ while the dimensionless parameters 
5/q and K represent the relative coupling constants of the planar quantum exchange interaction that 
respect/break the U(l) pseudospin symmetry and of another interaction that breaks the 27r -rotation 
symmetry of the pseudospins about the (111) direction. Figure 2 b shows the magnetic neutron- 
scattering cross-sections calculated at 0.3 K within the random phase approximation (RPAj^. A 
set of four adjustable parameters (J, S, q, K) = (0.68 K, —0.8, 0.2, —1.0) reproduces the experi- 
mental profiles shown in Fig. 2 a. (Any change of 6, q, or K by ±0.1 gives worse results.) The 
signs and ratios of the coupling constants obtained from the fitting are comparable to those ob- 
tained microscopically with a reasonable ratio Vpfj,/Vpfa- = 0.2 of the Slater-Koster parameter^^. 
The values of {S, q, K) are also nearly within the error bars of the recent spin-wave analysis of 
high-field inelastic neutron-scattering resultd^ (Supplementary Methods). The positive coupling 
constant J for the Ising interaction favours the 2-in, 2-out spin-ice manifold. (Note that the am- 
plitude of J could be underestimated here because of an overestimation of the onset temperature 
for the instability in the RPA.) Accordingly, a macroscopic number of states are nearly degen- 
erate, in spite of the additional interactions. In the theoretical calculations, the intensity and the 
sharpness of a remnant of the pinch-point singularity increase with decreasing temperature in the 



high-temperature paramagnetic phase (Supplementary Figure S4), in quaUtative agreement with 
our experimental results. These observations support the scenario that the high-temperature phase 
of Yb2Ti207 in the temperature range Tc < T < 2 J is well described as a classical magnetic 
Coulomb liquid which hosts the pinch-point features and monopolar spinons that are deconfined 
with a finite excitation energy. 

First-order transition to a ferromagnetic state In Fig. 2 e we present the temperature depen- 
dence in the sum of the SF and NSF magnetic neutron-scattering cross-sections across the [111] rod 
through (1.5, 1.5, 1.5). The broad intensity increases gradually on cooling from 0.7 to 0.3 K, but 
falls abruptly at 0.2 K, signalling a sudden change in the spin correlations. Actually, the neutron- 
scattering profile at 0.04 K (Supplementary Figure S3 d) clearly shows that the [111] -rod scattering 
loses intensity. We have measured the (111) Bragg peak intensity as a function of temperature on 
warming and cooling, with the incident neutron spins polarised along the X direction. Figure 3 a 
shows the dependence of the sum of the SF and NSF channels for small angular deviations from 
the wavevector (111) (rocking-curve scans) both above and below Tc with the data collected on 
warming the sample. By subtracting the 1 K data, the magnetic components can be clearly ob- 
served below Tc (Fig. 3 b). The integrated intensities for the SF and NSF channels (Fig. 3 c) and 
for the sum of the two (Fig. 3 d) show a remarkable hysteresis around Tc, indicating a first-order 
magnetic phase transition. It is possible that the magnetic scattering intensity along the diffusive 
[111] rod is removed to form the magnetic Bragg peaks below Tc. Well above Tc the flipping 
ratio, (the ratio of the NSF to SF cross-sections), at (111) gradually decreases on cooling and then 
falls more rapidly to unity around Tc (Fig. 3 e). This indicates that the incident neutron spins are 
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already partially depolarised in the high-temperature phase because of short-range ferromagnetic 
correlations, and are fully depolarised well below Tq, pointing to the emergence of macroscopic 
ferromagnetic domains. Furthermore, the first-order nature of the transition has also been con- 
firmed by a sharp anomaly in the specific heat of the sample, (Supplementary Figure SI), which 
closely resembles the earlier resulpi'. 

The emergent ferromagnetic order without translational symmetry breaking is also under- 
stood theoretically by the mean-field analysis (Supplementary Methods) of the same model at 
T = with the values of the dimensionless coupling constants (5, g, K) determined above. It is 
found that the mean-field ground state has a non-coplanar pseudospin order with the Ising compo- 
nent {S^) ~ 0.43 and the planar component ^Jl^sf)^^~+~{S^ ~ 0.25, as shown in Fig. 1 b. This 
is translated to a non-coplanar ferromagnetic order with a slight tilting from the [100] direction 
by ~ 1° and a moment amplitude 1.3yUB (Fig- 1 c), in fair agreement with the previous neutron 
diffraction results on the same single crystaP^. In this low-temperature phase, no soft spin excita- 
tions exist because of the absence of a continuous symmetry in the original model. This slightly 
noncoplanar ferromagnetic ground state does not alter with (5, g, K) within the accuracy of our 
fitting procedure. 

Discussion 

From the viewpoint of gauge theor)^"^, the high-temperature phase of our system is effectively 
described as a Coulomb phase, where bosonic fractionalised quasi-particles, "spinous", carrying 
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magnetic monopole charges are deconfined without a static non- uniform proliferation of monopole 
charges and host fluctuating U(l) gauge fields responsible for the analogous "quantum electrody- 
namics". These magnetic monopolar spinous can undergo a Bose-Einstein condensatiorP' to form 
the classical ferromagnetic momenP^ (Fig. 1 a). This pins the global U(l) phase of the spinon 
fields, and breaks the nontrivial emergent U(l) gauge structure. This condensation of matter 
(spinon) fields coupled to gauge fields occurs through the Higgs mechanism^. Thus, the ferro- 
magnetic state can be viewed as a Higgs phase of magnetic monopoles. 

Our observations support a novel picture that magnetic monopole charges are carried by frac- 
tionalised bosonic quasi-particles in the high-temperature magnetic Coulomb liquid. At present, 
the disappearance of the low-energy gapless excitation spectrum, unique to the Higgs phase, is 
only suggested by the long magnetic relaxation times seen in this and earlier studied, and must 
be confirmed by future experiments, which could also identify the spectrum of gapped spin wave 
excitations playing the role of Higgs bosons. It will also be fascinating to search for the quantum 
phase transition between the magnetic Coulomb phase and the Higgs phase with proper control of 
the material by pressure or chemical substitution. 

Methods 

Single-crystal growth Single crystals of Yb2Ti207, which were also used in the previous pub- 
lished workP^, were prepared by the floating zone method. Stoichiometric quantities of Yb203 and 
Ti02 powder were mixed, pressed into rods, and sintered at 1150°C for 24 hours. Using these 
rods, single crystals were grown in air at a rate of 1.5 mm/h. The crystals had a typical diameter of 
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~ 6 mm and a length of 20 mm. Powder X-ray diffraction measurements on a pulverised part of 
the single crystal showed no appreciable amount of any impurity phase. 

For comparison with other experimental work which reported the absence of long-range 
magnetic order, the different methods of sample preparation are worth mentioning. Only two re- 
ports are available in the literature. According to ref.|251 Gardner et al. synthesised poly crystalline 
samples of Yb2Ti207 by firing stoichiometric amounts of Yb203 and Ti02 at 1350°C for several 
days, with the quality checked by X-ray diffraction. According to ref. dH Ross et al. synthesised 
single crystals of Yb2Ti207 using the same floating zone method as our work, but in 4 bars of 
oxygen and at a faster growth rate of 5 mm/h. Single crystals prepared by Yaouanc et al. were 
also grown by the floating zone method but at an even faster growth rate of 8 mm/h in air with and 
without a subsequent heat treatment for 24 hours at 1 100°C under oxygen flov^'^. 

Characterisation of different single-crystal samples Here, we discuss the controversy over the 
presence or absence of a first-order phase transition in Yb2Ti207. We then demonstrate from 
measurements of the specific heat and the extended X-ray absorption fine structure (EXAFS) that 
the quality of the samples significantly alters the low-temperature behaviour of this material. In 
particular, we show that the long-range ferromagnetic order realised through the first-order phase 
transition in the single-crystal sample reported in the main text is easily removed by reducing the 
Yb content within the sample. 

First, let us briefly summarise the previous experimental reports that catalogue the strong 
sample dependence of the physical properties of Yb2Ti207. The Curie- Weiss temperature, 0cw> 



11 



varies considerably from sample to sample. 6cw ~ 0.53 K for our sample reported on in the 
main text and in ref. '24', 0.40 K in ref. |2ll 0.59 K in ref. 36, and 0.75 K in ref. |22l A significant 
sample dependence is also observed in the low-temperature specific heat in both powder and single- 
crystal Yb2Ti207 samples EHIIII, in addition, according to Gardner et al. a magnetic Bragg peak 
was observed at (1 1 1) with polarised neutron-scattering measurements on polycrystalline samples, 
although neutron spin depolarisation was not observed and static ferromagnetic order was ruled 
out. 

In an attempt to rationalise these conflicting reports, we have studied three single-crystal 
samples of Yb2Ti207 using specific-heat and EXAFS measurements. Sample A was prepared us- 
ing the method described above and was the sample used for both our polarised neutron- scattering 
measurements and the previous unpolarised measurements'^. Sample B was produced from a rod 
sintered at 1350° rather than 1 150° and was grown in air using a faster growth rate of 5 mm/h (cf. 
ref. [m who also used a faster growth rate but grew in 4 bars oxygen). Sample C was grown from 
another rod which was prepared in the same manner as sample A. However, the molten zone was 
less stable during the crystal growth. 

Supplementary Figure SI shows the temperature variation of the specific heat for these three 
samples. Sample A has a sharp singularity at 0.214 K. This is consistent with the observation 
made in our neutron- scattering experiments of the first-order phase transition to a ferromagnetic 
state. Similar sharp anomalies were also observed for polycrystalline samples in both the early and 
recent measurements EnSl. However, neutron- scattering experiments have not been reported for 
powder samples showing such an anomaly. For sample B, the specific-heat anomaly is broadened 
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significantly, signalling the disappearance of the first-order phase transition. The temperature de- 
pendence of the specific heat of sample B closely resembles those of the single crystals used 
for unpolarised neutron scattering experiments l ^^'^^'^^ i In sample C, the peak in the heat capacity 
disappears completely. 

The crystal structure of these three samples has been examined by EXAFS at the Yb L^-edge. 
Supplementary Figure S2 shows the radial distribution functions of these samples extracted from 
the Fourier transform of EXAFS oscillations, shown in the upper inset. The oscillations become 
increasingly damped for sample A (black), B (red) and C (green) respectively, indicating either a 
higher level of disorder or the appearance of more vacancies for sample A through to sample C. The 
Fourier transformed radial distribution functions favours the vacancy scenario. Near-edge spectra 
reveal only the Yb^+ charge state (see the middle inset). More oxygen vacancies, particularly at the 
0(2) site located at the centres of Yb40 tetrahedra, are created than Yb and Ti vacancies in samples 
B and C. These observations suggest vacancies of all three elements are created cooperatively to 
maintain charge neutrality. Comparing the EXAFS data with the specific heat curves, one can 
conclude that the higher the quality of the sample, i.e. the closer the sample is to stoichiometry, 
the sharper the magnetic transition. In particular, the low-temperature ferromagnetism discussed 
in the main text is intrinsic to the most ideal single crystals such as sample A. 

Polarised neutron-scattering experiments Measurements were performed at the high-flux po- 
larised diffuse neutron-scattering spectrometer DNS, FRM-II (Garching, Germany). A ^He/^He 
dilution refrigerator insert with an Oxford Instruments cryostat was used for the experimental tem- 
peratures from 0.04 to 1 K. A neutron wavelength of 4.74 A was chosen for all the experiments. 
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The [1-10] direction of the crystals was aligned perpendicular to the horizontal scattering plane so 
that the (hhl) reciprocal plane can be mapped out by rotating the sample. The neutron polarisation 
at the sample position was aligned along the [1-10] direction of the sample, i.e. the Z-direction of 
the chosen experimental coordinate system (Z-direction polarised neutron scattering). Within this 
setup, the SF and NSF scattering cross-sections are given by 



respectively, where M^yM^y and M^^M^^ are the components of the magnetic scattering cross- 
section in and out of the (hhl) scattering plane respectively, with Y being perpendicular to the 
scattering wavevector Q. Isi is the total nuclear spin incoherent scattering cross-section, and 
A^*A^ is the nuclear coherent scattering cross-section. In general, the magnetic scattering cross- 
sections MJ^yM^y and M]_zM±z are not identical. Therefore, they can provide information on 
the anisotropy of the magnetic correlations present in the system. In order to confirm the observed 
magnetic scattering, X-direction polarised neutron scattering has also been carried out. In this 
polarised neutron- scattering setup, the neutron polarisation is parallel to the scattering vector Q, 
and can be used to separate the magnetic scattering contribution to be mapped in the SF channel, 
and nuclear contribution in the NSF channel; 




(4) 



(3) 




(5) 



(6) 



For details, see refs [39ll40l 
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Figure 1 : Schematic phase diagram as a function of temperature T and the relative strength 
6 of the U(l)-symmetric planar exchange interaction and the Ising exchange, and the hypo- 
thetical low-temperature ordered structures of Yb2Ti207. a, Schematic phase diagram show- 
ing the first-order Higgs transition between a Coulomb liquid phase and a Higgs phase of magnetic 
monopoles. The other two model parameters q and K in Eq. ^ are assumed to be negligibly small 
for dipolar spin ice (Dy/Ho)2Ti207, while they are finite for Yb2Ti207 as obtained in the present 
work. Monopoles (blue balls) and anti-monopoles (red balls) are illustrated for both the phases. 
In the magnetic Coulomb liquid phase (yellow), magnetic monopoles are carried by pseudospin- 
1/2 fractionalised gapped spinon excitations out of quasi-degenerate spin-ice manifold, obeying 
a Coulombic law. In the Higgs phase (cyan), monopolar spinous are condensed to form local 
magnetic dipole moments (arrows) showing ferromagnetic long-range order, b, c. Hypothetical 
ferromagnetically ordered structures of the pseudospins (b) and the magnetic moments (c) in the 
low-temperature Higgs phase. The finite planar components of the pseudospins are ascribed to a 
condensation of monopolar spinous in the U(l) gauge theorjP'^. 
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Figure 2: Magnetic diffuse neutron-scattering profiles in the (hhl) plane with the incident 
neutron spins polarised along the Z direction, a. Experimentally observed profiles for the SF 
(left), NSF (middle) and total (right) magnetic diffuse neutron scattering cross-sections above Tq 
at T = 0.3 K. The background levels for the diffuse scattering have been subtracted by making 
use of the data measured with the incident neutron beam polarised along the X direction paral- 
lel to the scattering wavevector q (see Supplementary Methods and Supplementary Figure S3). 
White circles mark the Bragg spots where the background subtraction scheme is less reliable. The 
measurements were performed in the energy-integrated diffraction mode up to an energy of en- 
ergy 3.64 meV. b. Theoretically calculated magnetic profiles for SF (left), NSF (middle) and total 
(right) magnetic diffuse neutron scattering cross-sections at T = 0.3 K. The profiles were obtained 
as the equal-time correlators in the random phase approximation (Supplementary Methods) to the 
pseudospin-1/2 model, Hd + Hge- The form factor is not included, c and d. Neutron-scattering 
intensity for the SF channel along [11 — 1] and [111] cuts, respectively, across the remnant of the 
pinch point (111) atT = 0.3 and 1.5 K. The data were measured with a higher resolution in re- 
ciprocal space than the full (hhl) profiles. One or two points marking a high level of intensity at 
and close to the (111) point correspond to the nuclear Bragg peak. The lines are a guide to the eye. 
e. Temperature dependence of the total cross-section along [11 — 1] cut through (1.5, 1.5, 1.5), as 
indicated in the inset. 
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Figure 3: Hysteresis in the Bragg peak intensity and flipping ratio at (111) on cooling and 
warming across Tc, measured in the X-direction polarised neutron-scattering experiments. 

In order to ensure that the system was in thermal equilibrium, measurements well below Tc were 
performed after waiting longer than the magnetic relaxation time, which is of the order of 2 hours 
at 0.03 K (ref.l24l). a. The rocking-curve scans of the (111) Bragg peak for the sum of the SF and 
NSF channels measured while warming the sample. The horizontal axis represents the deviation of 
the wavevector from (111) with their amplitude being unchanged. The lines are a guide to the eye. 
All the peaks are instrument resolution limited, b. The magnetic contribution to the (111) Bragg 
peak is obtained by subtracting the data measured at 1 K, which is well above Tc, in the warming 
sequence, c. The growth of the integrated intensity for the (111) Bragg peak begins abruptly below 
Tc with a discernible hysteresis in both the SF and the NSF channels, d, The sum of the SF and 
the NSF (111) Bragg peak intensities versus temperature showing a hysteresis with temperature. 
The sum increases by almost 6% below Tc in agreement with the previous measurement^, e. The 
flipping ratio of the neutron spins at the (111) Bragg peak as a function of temperature. This ratio 
falls off steeply to unity below Tc with a clear hysteresis. The data were collected by first warming 
and then cooling. This sequence is indicated by the arrows in c, d and e. 
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Supplementary Figure 1 
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Supplementary Figure SI: Low-temperature specific heat versus temperature for three 
single crystals of Yb2Ti207. Sample A was used for the polarised neutron- scattering experi- 
ments reported in the main text. The heat capacity measurements were carried out in a Quantum 
Design Physical Property Measurement System (PPMS) using a pulse-relaxation calorimetry 
technique. 
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Supplementary Figure 2 
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Supplementary Figure S2: EXAFS radial distribution functions at the Yb La-edge from 
the Fourier transform of raw spectra in the upper inset. The near edge x-ray absorption 
spectra (XAS), shown in the middle inset, reveal only the Yb^+ peak at 8947 eV. The data are 
presented for crystals A (black), B (red) and C (green). 
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Supplementary Figure 3 
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Supplementary Figure S3: Raw polarised neutron-scattering cross-sections taken in the 
energy-integrated diffraction mode at DNS. a, b, and c, The SF, NSF, and total cross-sections, 
respectively, in the Z-direction polarised neutron-scattering measurements at T = 0.3 K. d. The 
total cross-section in the Z-direction polarised neutron-scattering measurements at To = 0.04 K. 
e and f. The SF and NSF cross-sections, respectively, in the X-direction polarised neutron- 
scattering measurements at Tx = 0.7 K. Sharp peaks at reciprocal lattice vectors (111), (002), 
(220) and so on are dominated by the nuclear contributions. The weak peaks at (1.5, 1.5, 0.5) 
and (112) are attributed to contamination from higher order scattering and are not intrinsic. 
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Supplementary Figure 4 




Supplementary Figure S4: Temperature dependence of the theoretical calculated equal- 
time correlator corresponding to the energy-integrated Z-polarised neutron-scattering 
cross-section in the SF channel, a, The intensity along the [11 — 1] cut through the (111) 
point, b, The intensity along the [111] cut through the (111) point. 
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Supplementary Methods 



Background subtraction 




Here, we explain how we derive the Z-direction polarised diffuse magnetic scattering profiles at 
0.3 K shown in Figure 2 a of the main text from the raw data shown in Supplementary Figure [S3| 
We estimate and subtract the background levels for the SF and NSF channels as follows. 

Because of the formation of a long-range ferromagnetic order and the loss of the [111] rod 
diffuse scattering intensity well below Tq K, the sum of the SF and NSF channels at To = 0.04 K 
can be roughly regarded as a temperature-independent background due to the nuclear coherent 
and incoherent scattering cross-sections; 

Total bg / rl \ 

=^-^+^-ni),„ , (SI) 

except around the Bragg peak positions where we suffer from a slight overestimate due to the 
magnetic Bragg-peak intensities. On the other hand, the SF and NSF cross-section profiles 
obtained at 0.04 K cannot be used as background levels, because the incident neutron spins are 
fully depolarised at 0.04 K, while they are only partially depolarised at 0.3 K, as is clear from 
Figure 3 e of the main text. 

Note that the spin flipping ratio is defined as the ratio of the NSF to SF scattering intensity 
at a Bragg-peak position (111), where nuclear coherent scattering cross-section N*N, which 
appears only in the NSF channel according to Eqs. (5) and (6) as well as Eqs. (3) and (4), 
dominates over the other contributions. Then, it is reasonable to assume that the SF and NSF 
cross-sections observed at a temperature T with the Z- and X-direction polarised neutron- 
scattering measurements take the forms 



Z-SF 



Z-NSF 




[MlyM^y)^ + Usi + r^' \{MIzM^z)t + N*N + \lsi 



1 + r^' 

{MlyM^Y)T + Ihi] + {MIzM^z)t + N*N + \lsi 



1 + r^' 

{MlyM^Y + MIzM^z)t + Ihi + \n*N + \lsi 



{MlyM^Y + MIzM^z)t + Ihi] + N*N + \lsi 



l + rn 



(S2) 



, (S3) 



(S4) 



(S5) 



Using Eqs. (|S4l ) and (|S5] ). the magnetic cross-section can be eliminated and the temperature- 
independent purely X-NSF background cross-section is given by 



at a temperature Tx- For a higher accuracy, it is better to take the data with a larger spin flipping 
ratio ttx^ and thus we take Tx = 0.7 K where ~ 14.97 (Figure 3 e of the main text). The 
raw data for the X-direction polarised elastic diffuse neutron- scattering experiments are shown 
in Supplementary Figures [S3] e and|S3]f. In fact, the q-space region where the X-direction 
polarised experiments have been performed does not cover all the reciprocal space required for 
this analysis. Nevertheless, the profile for Eq. (|S6l) obtained from Supplementary Figures [S3] e 
andlS3]f at = 0.7 K is reproduced to an accuracy of 10%, except at the Bragg-peak positions, 

/ , \ total I I 

by multiplying the profile for y^j^ at Tq = 0.04 K shown in Supplementary Figure[S3]d by 
an overall scale factor 0.6. This allows us to simulate the NSF background profile for Eq. (lS6]) 
in the whole required area of the [hhl) plane. 

Now, we can express the background cross-sections for the SF and NSF channels in the 
Z-direction polarised neutron-scattering experiments as 




Z-SF 



bg []v*X + Isi] - (1 - r^') \n*N + \lsi 



l + r-r, 



(S7) 
(S8) 



where [N*N + Igi] and 
flipping ratio at T = 0.3 



N*N + llsi 



are given by Eqs. (iSTl) and ((S6]) . respectively. The spin 
C is given by = 8.7 (Figure 3 e of the main text). 
In principle, we should also be able to eliminate a small (11.5%) mixing of MJ^^M^^ in the 
SF channel and of M^y^±y in the NSF channel at 0.3 K. However, this effect is found to be 
small compared with the accuracy of the background subtraction. 



Mean-field theory and random phase approximations 

The classical mean-field (MF) Hamiltonian is given by 



Hmf - ^ J2 ('^^-q,i)K'^i,j('^^q,j) 



q fj.,i/=x,y,zij=0 



(S9) 



where h'^'^^ j is obtained from the Fourier transform of the full Hamiltonian {Hd + Hge) given in 
the main text^^'^^, and i and j denote indices for the Yb sites at i? + aj (i = 0, ■ ■ ■ , 3) with fee 
lattice vectors R and ao = (0, 0, 0), ai = a(0, |, ^), a2 = 0, and as = a(^, |, 0), and 
the cubic lattice constant a. The local coordinate frames are taken as 



S^-^(U.-2) 
~ (1,-1,2), e? 



-,y 



-^(-1,1,0), e^ = -L(l,l,l^ 



-^(-1,-1,0) 



^(1,-1,-1), 



(SIO) 
(Sll) 
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e2 = -;^(-l,l,2), el = --J= (1,1,0), = -L(-l, l, -l), (S12) 

e? = -^(-l,-l,-2), e| = --^ (1,-1,0), e^ = -L(-l,-l,l), (S13) 

to maximally exploit the symmetry of the system. In particular, the z axis points to the (111) 
direction of the Ising magnetic moment. This in turn introduces the following simple bond- 
dependent phases (^^y = 4>ij = 4>ji in the two superexchange terms proportional to q and K 
which break the planar U(l) pseudospin symmetry; 

27r 27r 

001 = 023 = g-, 002 = 031 = -g-, 003 = 012 = 0, (S14) 

where (i, j) is a pair of sublattice indices for r and r' . 

The Ewald summation technique was employed to calculate the long-range magnetic dipolar 
interaction. Then, the eigenvalue problem is solved for each q as 

T.KljK,n,-^.,nK,n,. (S15) 

V 

where n labels the twelve eigenvalues/eigenstates. The magnetic neutron-scattering cross- 
sections at a temperature T are then calculated with the random phase approximation (RPA)^^'^^ 
as 

n i j=l + ^ 



where 



fZ,^ ^^,\Z.{q^ {g,et x q))] [Z • (q x {g^e) x q))] , (S17) 
with Z = (1, -1, 0)/\/2 for the NSF channel and 

for the total. The SF channel is obtained as the difference between the above two. 

In our case of J > 0, while 5 > leads to profiles analogous to the nearest-neighbour 
Ising exchange or dipolar spin ice^", a choice of 5 < produces various profiles depending 
on the other relative coupling constants q and K. The best fit to our experimentally obtained 
and analysed magnetic scattering profiles at 0.3 K, shown in Figure 2 a of the main text, was 
obtained with (J, 5, g, K) — (0.68 K, —0.8, 0.2, —1.0), as shown in Figure 2 b of the main text. 
This set of parameter values is compared to the result from the spin-wave fitting^", which is 
translated into J = 2.0 ± 0.5 K, 5 = -0.6 ± 0.3, g = 0.6 ± 0.3 and K = -1.6± ~ 0.4 in 
the present notation. Our exchange energy scale is smaller by a factor of 2-3. The discrepancy 
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in J may be partly attributed to an overestimate of the temperature at which an RPA instability 
occurs, ~ 0.29 K. This gives 6'cw ~ 0.29 K. Scaling the energy scale by a factor of 2 reproduces 
the experimental observation 6'cw ~ 0.53 K (ref. 24). 5 and q are nearly within error bars of the 
analysis given in ref. 20. However, we note that if we take \5\ to be comparable to q and/or take 
a much larger value of \K\, e.g., K = —1.2, following ref. 20, the intensity at the scattering 
branch from (1 1 1) to (220) in the SF channel is significantly suppressed, in disagreement with 
our experiments. The nearest-neighbour coupling constants from Hjj are of the order of 0.01 K, 
namely, one or two orders of magnitude smaller than the superexchange couplings of our choice, 
though they give rise to ~ —0.2 K in the energy eigenvalue of the mean-field Hamiltonian. This 
might cause other minor differences. As also noted in ref. 20, the parameters 5, q, and K 
we have obtained are quite different from the estimates by Thompson et al^ who analysed the 
Hamiltonian including all the crystal field states but only the same number of coupling constants 
without any multipolar interactions. 

In Supplementary Figure|S4]a and b, we show the temperature dependence of the calculated 
neutron- scattering intensity in the SF channel (without taking account of the form factor) along 
the [11 — 1] and [111] cuts through the (111) point, respectively. It is clear that the intensity 
along the [11 — 1] cut is not only enhanced but also sharpened with decreasing temperature 
down to 0.3 K, while that along the [111] cut remains almost featureless even at 0.3 K. This 
is consistent with the picture that the high-temperature phase of this model is described as a 
Coulomb phase characterised by a remnant of a pinch-point singularity. 

The mean-field ground state was obtained from the same scheme as employed above. We 
searched for the lowest-energy states among the classical energies per site, 

Ecl,q,n MaXiY.t,\(t)q,n,i?' ''^^^'^ 

Then, the MF ground state is found by approaching the (000) point from the [100] direction, 
resulting in the noncoplanar pseudospin structure and the nearly coUinear ferromagnetic struc- 
tures shown in Figures 1 b and c of the main text, respectively. 

Note that the possibility of an ordered spin-ice state^"'"^^'"^^ is ruled out, since it should create 
magnetic Bragg peaks at (100), in contradiction to our findings and cannot explain the neutron 
depolarisation due to the macroscopic magnetic moment. 
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